Generation of electric currents in the 
chromosphere via neutral-ion drag 
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ABSTRACT 

We consider the generation of electric currents in the solar chromosphere where the ionization 
level is typically low. We show that ambient electrons become magnetized even for weak magnetic 
fields (30 G); that is, their gyrofrequency becomes larger than the collision frequency while ion 
motions continue to be dominated by ion-neutral collisions. Under such conditions, ions are 
dragged by neutrals, and the magnetic field acts as if it is frozen-in to the dynamics of the 
neutral gas. However, magnetized electrons drift under the action of the electric and magnetic 
fields induced in the reference frame of ions moving with the neutral gas. We find that this 
relative motion of electrons and ions results in the generation of quite intense electric currents. 
The dissipation of these currents leads to resistive electron heating and efficient gas ionization. 
Ionization by electron-neutral impact does not alter the dynamics of the heavy particles; thus, 
the gas turbulent motions continue even when the plasma becomes fully ionized, and resistive 
dissipation continues to heat electrons and ions. This heating process is so efficient that it can 
result in typical temperature increases with altitude as large as 0.1-0.3 eV/km. We conclude that 
this process can play a major role in the heating of the chromosphere and corona 
Introduction equally difficult to understand. 



The detailed physical mechanism of coronal 
heating is not yet well understood. A number 
of fundamental questions remain that challenge 
theoretical descriptions and the interpretation of 



IParkerl (|l988h proposed the idea that the solar 



corona could be heated by the episodic dissipa- 
tion of energy at many small-scale tangential dis- 
continuities arising spontaneously in the coronal 

magnetic field, as it becomes braided and twisted 
observational data (see, e.g. Khmchuk, 2006, or , , , i_ • <? , • m 

q j onno t ■ \ t j by random photospheric footpomt motions. Ihese 



Walsh & Ireland, 2003, for recent reviews). In ad- 
dition, the heating of the chromosphere require; 
much more energy, and this related process is 



. events - sudden changes of the magnetic field 
dition, the heating of the chromosphere requires , , ,, ,. ,, . , , 

m „ ra a 4-u- i „ • topology - hypothetically result m plasma heat- 

ing and the acceleration of non-thermal particles. 
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ration for this wa s the discovery of the hard X- 
ray microflares by Lin et all ( 1984 ) ; the energy of 
one nanoflare was to be roughly 10 -9 times the 
energy of a major flare, and thus orders of mag- 
nitude weaker than even the microflares. Parker's 
idea stimulated intensive searches for any obser- 
vational signatures o f nanoflares ( Cargi I Il994t 
Sakamoto et al.l 120081 IVeksteinl 120091 ) and their 
possible contributio n to the overall energy budget 
of the solar corona (|Klimchukll2006( ) . 

Microflares were first detect ed in hard X-ray s 

" f|l984h . 



Lin et al 



in a balloon experiment by 
The subsequent development of new instrumenta- 
tion allowed the multi-wave satellite and ground 
based high-resolution observations of smaller-scale 
(about a thousand of kms or even smaller) and 
lower energy phenomena. Soft X-ray imaging 
revealed abundant m icroflares in active regions 
(ShimizuetaL and RHESSI observations 

found that virtually all of a sample of some 25,000 
hard X-ray microflar es occurred in active regions. 
Krucker et al.l ([19971 ) found flare-like brightenings 
in areas of the quiet Sun, and observations at 
EUV wavelengths (e.g., the "blinkers" of Har- 
rison, 1997) reveal bursting activity above the 
magnetic network borders. Similar phenomena 
that form small X-ray jets at the limb were re- 
ported by Koutchmy et al. (1997). From these and 
other EUV observatio ns (jBerghmans et al" 1998t 
Benz fc Kruckerl ll998). if not the hard X-rays, we 
conclude that Parker's idea of episodic heating of 
the apparently steady quiet corona should not be 
discarded, even though n o convincing e vidence for 
the required steepening (|Hudsonlll99ll ) of the en- 
ergy distribution function has yet been presented. 

However, the idea of coronal heating via tan- 
gential discontinuities that arise spontaneously 
(nanoflares) does not address two important ques- 
tions. Namely, where does the excess magnetic en- 
ergy come from, and what was its original source? 
Parker's analysis considered ideal MHD statisti- 
cal eq uilibria that contain m ultiple discontinuities. 
Later iRappazzo et al. I (|2008h showed that a large- 
scale MHD energy source perturbed by slow mo- 
tions on its boundary, supposed to be induced 
from the photosphere, results in the generation of 
a Poynting flux. This drives an anisotropic turbu- 
lent cascade dominated by magnetic energy. The 
result looks similar to Parker's tangling of mag- 
netic field lines but the small-scale current sheets 



(which replace the tangential discontinuities) are 
continuously formed and dissipated. In this mod- 
ification of the initial scenario the current sheets 
are the result of the turbulent cascade. The ini- 
tial energy reservoir in such a view is contained in 
large scale magnetic-field structures. 

Here, we discuss another possibility: the direct 
generation of relatively small-scale electric cur- 
rents by neutral gas motions. Clearly, a huge 
energy reservoir exists in the form of turbulent 
motion of neutral gas at and beneath the pho- 
tosphere, supported by the underlying convection 
zone. It is widely accepted that this energy can 
be partially transformed into the excess magnetic 
energy in the chromosphere and corona. How- 
ever, there is no quantitative model that describes 
the physical mechanism of such energy transfer. 
The development of a model that describes the en- 
ergy transfer from quasi-neutral gas motions to the 
magnetic field in a step-by-step manner inevitably 
raises questions about the spatial and temporal 
scales at which such a transfer can occur and about 
its location. Recent observations provide strong 
indications that the energy reservoir is indeed the 
dynamic turbulent photosphere, and that the en- 
ergy transfer from the turbulent gas motions to 
various kinds of trapped and transient magnetic 
field oscillations takes place at the chromospheric 
level. 

Analysis of the high resolution spatial (^150 km 
at the Sun) and temporal (few sec) data obtained 
by the SOT (Solar Optical Telescope) aboard 
Hinode (De Pontieu et al., 2007) revealed that 
the chromosphere is dominated by a multitude of 
thin (~200 km wide) dynamic, jetlike "type II" 
spicules. They are ejected upwards with charac- 
teristic velocities of 20-150 km/s and reach heights 
of 2000-10000 km before disappearing from the 
chromospheric passband (in this case, the H line 
of Ca n) . The type II spicules have quite short life- 
times of 10-300 s (according to the authors, most 
of them last less than 100 s) and many of them 
undergo substantial transverse displacements of 
the order of 500-1000 km. Moreover, the large- 
scale long-living spicules display oscillatory mo- 
tions in the direction perpendicular to their own 
axes. Since the spicule structure can be taken to 
outline the direction of the magnetic field, this led 
the authors to the conclution that the observed 
motions indicate Alfvenic oscillations. 
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Furthermore, the spatio-temporal variations of 
the chromosphere have always revealed the great- 
est complexity, and prominences also consist of nu- 
merous threadlike feature s with strong an d mixed 
flows along these threads ( Lin et al]|2005l) . Obser- 
vations with SOT have provided an exceedingly 
variable and dy namic picture of thes e flows and 



field structures (|Okamoto et al.ll2009t ). The SOT 



chromospheric data are in the H-line (Ca n), show- 
ing plasma at roughly 2 x 10 4 K. The movie pre- 
sented by Okamoto et al (available on the Hinode 
web site http://solarb.msfc.nasa.gov/) shows 
ubiquitous continuous motions along the promi- 
nence thread lines. The oscillatory motions ob- 
served might be interpreted in terms of propagat- 
ing or standing Alfven waves on the magnetic field 
lines that presumably compose the prominence. 
The typical transverse spatial scale of threads was 
found to be of the order of 600 km, with a charac- 
teristic length of the order of several Mm. The 
characteristic temporal scale was found to vary 
from 100 to several hundred seconds. 

SOT observations near the limb led to the 
discovery of another small scale dynamic phe- 
nomenon in the chromosphere: tiny chromo- 
spheric "anemo ne jets," na med for the similar 



X-ray features (jShibata et al.l 120071 ) . They also 



resemble larger-scale features well known in Ha 
data and called surges or sprays, which often oc- 
cur near sunspots and in association with flares 
or other transient activity (e.g. Rust, 1968). All 
of these observations can be considered as indica- 
tions of the generation of small-scale perturbations 
and oscillations of the magnetic field in the chro- 
mosphere, where the degree of ionization is still 
relatively small. 

2. The coupling problem 

In spite of greatly improved observational data, 
the physics of chromospheric and coronal heating 
is not yet well understood. While the underlying 
mechanism must be associated with the magnetic 
field, the details of how efficiently the energy of 
convective motion is transformed into magnetic 
energy in the solar atmosphere is an open ques- 
tion. The importance of this process should not 
be understated. It is relevant not only to coro- 
nal heating and the acceleration of the solar wind, 
but also to the formation of the initial spectrum 



of wave turbulence introduced into the solar wind. 

We propose a mechanism to generate electric 
currents that can be viewed as Alfven or magne- 
tosonic waves that result from strong ion-neutral 
drag. We assume that photospheric motions 
are turbulent, and consist of both compressional 
and rotational flows with energies exceeding that 
which is necessary to drive the observed oscilla- 
tions of suspended threads or jets in the chro- 
mosphere. The question then becomes how these 
oscillations can be transported from the photo- 
sphere to the chromosphere, and how the flow of 
neutrals is transferred into the motion of charged 
particles and, ultimately, into the generation of 
electric current. 

There are two physical processes that are nec- 
essary to convert the energy of the neutral gas 
motion into the magnetic field oscillations. Pho- 
tospheric motions of charged particles are domi- 
nated by frequent electron-neutral and ion-neutral 
collisions, so that ions and electrons tend to follow 
the neutral gas motion with the zero net electric 
current. 

Since density falls off rapidly with height in 
the solar chromosphere, the electron collision fre- 
quency with neutrals and ions decreases to the 
point where it becomes smaller than the elec- 
tron gyrofrequency. However, the ion- neutral col- 
lision frequency at this height still remains sub- 
stantially larger than the ion gyrofrequency, since 
ion-neutral and electron-neutral collision rates are 
not substantially different (see, e.g., De Pontieu 
et al., 2001). Therefore, the motions of ions and 
electrons differ. Electrons tend to move along the 
magnetic field lines and drift due to E x B in 
the transverse direction, while the ions continue 
to move together with the neutrals. This differ- 
ence results in the generation of electric currents. 

To view this situation from a different perspec- 
tive, if ions move together with neutrals (due to 
the strong drag) with some angle to the back- 
ground magnetic field, an inductive electric field 
v x B appears in the reference frame of the 
plasma. This induced electric field inevitably gen- 
erates electric currents that can be calculated us- 
ing the plasma conductivity tensor. When elec- 
trons and ions are both demagnetized, this tensor 
reduces to the simple scalar conductivity, and the 
resulting current is very small. When electrons 
become magnetized, the motions of electrons and 
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ions become decoupled, and the efficiency of the 
current generation is substantially higher. 

Thus, we conclude that there is an efficient 
chromospheric dynamo operating in the layer be- 
tween the level of electron magnetization and the 
height where either the degree of ionization be- 
comes high (comparable with unity) or where the 
ions become magnetized. Forced oscillating cur- 
rents are generated most effectively when the in- 
ductive electric field has characteristic frequencies 
and wavelengths close to the eigenmodes of the 
system, i.e. the MHD waves in the magnetized 
weakly ionized plasma. At higher altitudes, where 
ions are also magnetized, the motion of ions and 
electrons can still differ. However, under these 
conditions, both components will mainly move 
along the magnetic field lines and perform drift 
motions across the field. Thus, the current is di- 
rected mainly along the magnetic field lines, while 
in the intermediate region it can be generated in 
an arbitrary direction. 

In summary, understanding how currents are 
generated by the turbulent motions of neutral gas 
in the photosphere can be broken down into two 
steps. First, we must understand how the energy 
and vorticity of photospheric motions is trans- 
ported upward to the level where electron and ion 
motions become decoupled. Second, we must de- 
rive a self-consistent system of equations that de- 
scribe the inductive electric field and the result- 
ing deformation of the background magnetic field. 
Once we obtain the electric field and estimates of 
current density, we can then determine the effi- 
ciency of electron heating due to collisions and can 
investigate the efficacy of other possible heating 
mechanisms. 

3. Atmospheric models 

We appeal to standard semi-empirical model at- 
mospheres for approximate values of the physical 
parameters of the plasma in the pho tosphere and 
chromosphere. iFontenla et al. (2009) provide a re- 
cent series of eight such models representing fea- 
tures such as the quiet Sun, faculae, and sunspots. 
Such models are based on macroscopic radiative 
transport theory and are adjusted to recreate so- 
lar spectroscopic observations. They do not repre- 
sent the dynamics, nor the plasma physics, since 
computer technology thus far allows at best only 



an M HD approach to the physics (e.g., ISteiner 
20071 ). Figure [1] compares co llision frequencies 



from the IFontenla et al.l ([2009) quiet-Sun model 
(Model 1001) with the gyrofrequencies of electrons 
and ions. For this model we assume a constant 
magnetic field o f 30 G (e.g. lHarvev et al. 2007 ; 
iLites et al.ll2008h . 




-0.5 0.0 0.5 1.0 1.5 2.0 2.5 
Height, Mm 



Fig. 1. — Collision frequencie s in the quiet-Sun 
model of Fontenla et al.l ( 20091 ). Blue (red) curves 
represent proton (electron) collision frequencies, 
with the lower lines showin g neutral rates only 
(e.| ' 



De Pontieu et a l. 2001); the dotted horizon 



tal lines show the gyrofrequencies (ions and elec- 
trons) for an assumed field strength of 30 G. 

The Fontenla models cover a range of physical 
features in the solar atmosphere, and were origi- 
nally intended for irradiance modeling. Here, we 
use these models as a guide to set a range of pa- 
rameters that are reasonably consistent - in this 
limited theoretical framework - with solar struc- 
ture. Figure [2] shows the location of the chromo- 
spheric dynamo layer for the full range of models. 
Electrons become magnetized, for all models, es- 
sentially throughout the solar atmosphere. Thus 
the chromospheric dynamo layer begins close to 
the photosphere and extends high into the chro- 
mosphere. This is true for all of the Fontenla 
models, ranging from sunspot umbra to bright fac- 
ula, for which the transition-region pressures range 
from 0.2 to about 2.0 dyne cm -2 . Note that in 
such models, substantial hydrogen ionization does 
not take place below the transition region between 
chromospheric and coronal temperatures. 

We emphasize that these models do not have 
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any plasma physics in them as such, and only serve 
as references at the order-of-magnitude scale. In- 
deed, the mechanism discussed in this paper would 
point to areas in which models of the solar atmo- 
sphere can be greatly improved. 

Dynamo layer 



Bright facula 
Sunspot penumbra 
Sunspot umbra 
Facula 

Plage (not facula) 
Enhanced network 
Quiet-Sun network lane 
Quiet-Sun inter-network 
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Fig. 2. — Chromospheric dynam o layers for the 
variou s solar features modeled by iFontenla et al 



(|2009h . The field strengths assumed are 30 G for 
the quiet Sun (the lower two models), and 1500 G 
for the others. 

4. Energy transfer from the photosphere 
to the chromosphere 

We first consider the physics of mass, momen- 
tum, and energy transport in the lower chro- 
mosphere. In order to formulate a mathemati- 
cal description of the problem we need to take 
into account physical conditions between the pho- 
tosphere and lower chromosphere, as described 
above. 

We begin with an isothermal atmosphere, i.e., 
a gas having a constant sound speed and strati- 
fied by a uniform gravitational force acting in the 
negative z-direction, for which 



po (z) = p ph exp {-z/H) . 



(1) 



Here p p h is the density at the level of the photo- 
sphere corresponding to z = 0, and H is the den- 
sity scale height. The gas motions are described 
by the equation of continuity 



dp 
dt 



+ V-{(po + p)^} = 0, 



where the density is presented as the sum of hy- 
drostatic equilibrium density and a perturbation 
(note that in hydrostatic equilibrium there are no 
flows) . 

By taking into account the equilibrium density 
dependence upon height this can be re-written as 



1 (£. 

dt \p 



dv x dv y 
dx dy 



d 



1 



dz H 



H Pq \Pq J 

where we have separated the linear and nonlinear 
parts and put the latter on the right-hand side. 
The equation of motion 



{P + Po) ( 



dlf 



dt 

takes the form 



(^V)^) = -VP-(p + p )G5 



— + V — -z- — + — 

dt po H po po 

= -( vV )v- £-& + (vV)v). (3) 
po dt 

We shall also use the equation of state taken in 
the simplified form of the equation for adiabatic 
motion, 

dP = i P ph d 
dt 



This is 



Pph dt 
d 



(P + Po) ■ 



dP_ 

~dt ~ a dt 
It also can be re-written for perturbations as 



d_,P_. 
dt [ p ' 



r 2 d ( p u 7 ~ 1 °s 



' dt po' 



7 H 



cSK^v)t-5(t-)]-C^v)(^)+§(?)- (4) 



Po H po' 



"po H po' 



Our objective here is the analysis of propaga- 
tion of the photospheric perturbations upwards 
into the chromosphere, and for this purpose we 
analyze the characteristics of linear perturbations. 
To obtain the linear wave mode dependencies, we 
determine the eigenmodes of the linearized system: 

d P „ _v 1 „ dp _ v 1 

JL+S7.Tf---v z = => —J- = -V--&+—V,, 
dt po H dt po H 

(5) 
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d~Tf P IP p 

+ v(-) - z-(-) + (^)Gz = 0, (6) 
at Po H pa po 

and 

|(^) + ClV.l^-i^ = 0. (7) 
at po 7 

We choose standard solutions of the form: 
P, p,v ~ exp i(u;i — /c^rr — /c y y — fc 2 ,z), 

and analyze the dispersion relations that deter- 
mine the dependencies of the frequency uj upon 
components of the wave vector k . But our prob- 
lem is formulated as a boundary-value problem 
where we assume the perturbations and their time 
dependences to be prescribed at the boundary 
z = 0. We want to describe the evolution of per- 
turbations with height z. In this case the very 
same equations should be used to find out the 
dependence of the component k z of the fc-vector 
upon the frequency u> and other components k x , 
k y . Substituting dependencies defined above into 
the linearized equations one can obtain the disper- 
sion equation in the form 

[c 2 - lK g}{^ - Cl[(e + K 2 ) + ^^1)]}- 

H 7 

lG C 2 (k 2 + K 2 )(K+^- 1 —±)=Q, 
H 7 

where C*| = -yGH, k 2 = fc 2 + k 2 yl K = k z - jj . 



i , / W 4_ w 2 C 2 (fc2 + _^ ) + (7 _ 1)fc2;? 2 



k z = ± 

2H 



This analysis demonstrates that a wide class of 
perturbations satisfying the condition 

c 4 - w 2 C 2 (fc 2 + i l I ) + ( 7 -l)fc 2 G 2 >0 

increase with height, and the characteristic growth 
rate is 

This phenomenon, the "effective growth" of 
the perturbations with altitude in a hydrostatic 
gas equilibrium with exponentially decreasing 



density is well known in the terrestrial atmo- 
sphere. It plays an important role in the prop- 
agation of infrasonic perturbations induced by 
explosions in the atmosphere (e.g. Klostermeyer, 
1969a, 1969b) . The same physical phenomena 
play important roles in the solar photosphere 
and chromosphere and have been the subject 
of many observational, theoretical, and simula- 
tion studies. At low frequencies acoustic-gravity 
waves (purely hydrody namic) play a role (e.g. 



MHD modes 
tions (e.g 



Straus et al.l 12008. 200 9J), and at h igh frequencies 



e.g.. lHollwesj|l986l) The simula- 



Stein fc Nordlundl Il998h can include 



essential physics such as optically-thick radiative 
transfer. 

The results of the more realistic simulations 
highlight important features of neutral gas motion. 
For example, vorticity is primarily generated due 
to baroclinic forces, and tends to concentrate in 
tube-like structures whose widths are comparable 
to the numerical resolution. In addition, plasma 
heating depends on the formation of supersonic 
flows and shock-like structures. If shocks form 
systematically, they could well be an important 
factor in the physics of coronal heatin g. However 
according to lStein fc Nordlundl (|1998l ). shocks are 
indeed observed at the edges of integranular lanes, 
but are a rare occurrence. At any one time, su- 
personic flow occurs in only 3-4% of the surface 
area. 

The role of acoustic gravity waves has also been 
intensively studied experimentally. Recently new 
experimental studies combining the SOT/NFI and 
SOT/SP instruments on Hinode and the Michel- 
son Doppler Imager (MDI) on SOHO, with the 
aid of 3D computer sim ulations, were performed 
by IStraus et al. ( 2009t ). The authors came to 
the conclusion that the gravity waves are the 
dominant phenomenon in the quiet middle/upper 
photosphere and that they transport sufficiently 
more mechanical energy than the high frequency 
(>5 mHz) acoustic waves. In addition they con- 
clude that the acoustic-gravity wave flux is 3- 
5 times larger than the u pper-limit estimate of 
Fossum fc Carlssonl ((2005) . These observations 



together with the numerical models allow us to 
consider acoustic gravity waves as one of the most 
important sources from which energy may be 
transformed to the magnetic field, and then to 
heat via the action of electric currents. 



6 



We note that that the absolute magnitude of 
the density perturbation actually decreases with 
height, but it drops more slowly than the back- 
ground density. This results in the growth of the 
relative density perturbation, written as p/po hi 
our notation. On the other hand the velocity per- 
turbations really grow exponentially. An impor- 
tant issue here is the characteristic vertical scale 
of this growth. Taking the gas temperature to be 
5000 K we hnd 



H = 



k B T 
MG 



140 km. 



This linear analysis includes only the effect of the 
growth of the velocity and relative density pertur- 
bations with the altitude, hence, the perturbation 
scales do not vary with altitude. Of course, there 
are a wide range of phenomena, such as vortical 
flows, that are not accounted for in this formalism. 

For example, the vortex radius can decrease 
with altitude. Indeed, assuming conservation of 
angular momentum the increase of the velocity 
should result in the shrinking of the transverse di- 
ameter of the vortex. This effect is relevant for 
Rossby vortices in a multi-layer atmosphere and 
it is in a good agreement with the observation of 
Stein and Nordlund (1998) that the vorticity con- 
centrates in small scale tubes. 

5. Magnetic perturbations produced by 
the turbulent motions of the neutral 

gas 

Turbulent convective fluid motions at the pho- 
tospheric level extend upwards into the lower chro- 
mosphere. The latter is a weakly ionized plasma 
with a typical temperature of 6000-7000 K, neu- 

3 , and 



tral hydrogen density nn ~ 10 — 10 
electron density n e ~ 10 10 — 10 11 



14 



Under 



these conditions the frequency of the ion-neutral 
collisions is as high as v in ps 10 7 — 10 9 s _1 which, as 
will be confirmed below, produces quite a strong 
drag effect resulting in the bulk velocity of ions, 
Vi closely matching that of the neutral gas, V n . 
Then the magnitude of the generated electric cur- 
rent ~~] = n e e(Vi — V e ) « n e e(V^ — V^) is deter- 
mined by the electron bulk velocity V e . The latter 
is governed by the equation of motion for the elec- 



trons 

me—jj- = -e$+^ e xB]}-{v en +v ei )m e (v e -V n ), 

(8) 

where v en w 5 x 10 8 — 10 10 s 1 is the frequency 
of electron-neutral collisions, the major source of 
the electron drag. Here v e i is electron-ion collision 
frequency. 

Ions are supposed to be dragged by neutrals, 
thus the velocity of electrons can be expressed as 

en 

The physical processes we describe occur in a 
parameter range corresponding to the transition 
in electron motions from unmagnetized, collision- 
dominated, to magnetized, when initially \u>Be\ < 
(v en -\- v e i), and then with the growth of altitude 
|wBe| becomes larger than (v en + v e i)- We consider 
hereafter slow motions uj << (V e n+fei)> \^Be\ sim- 
ilar to those we described in the previous para- 
graph. In this case one can neglect the electron 
inertia in the left-hand side of the equation [5J It 
can then be simplified to obtain: 



-e{^+[V^x^] [fx~^]}+{ven+Vei)m e — = 

ne ne 

(9) 

Taking the curl and replacing V x B via 
one can obtain the final equation 



^— = V x [v£ x ih - —V x -[V x it x it] 
at epo n 



(v en + v ei )m e VxB 
— V x . 



e 2 p 



where j is replaced by 



J 



= Ivxl 



Mo 



(10) 



(11) 



Equation [TU] together with the system of equa- 
tions [2 [3l [4] make up a closed system describing an 
interaction of neutral gas perturbations with the 
ionized electron component that can result in gen- 
eration of electric currents. We neglect hereafter 
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the influence of electron motions on neutral gas 
motions by assuming that the degree of ionization 
is low. 

In what follows it is assumed that the chromo- 
spheric magnetic field can be represented as 



where Bq is a background field, which for the sake 
of simplicity is assumed here to be just a uniform 
field, while b is a relatively small field deformation 
caused by a prescribed flow of neutral gas with 
velocity V n and frequency oj. 

Then the linearized (with respect to b) version 
of Eq. [TU] takes the form 



t+idl ] -^(li-VWxt-d 2 Jl+i^)V 2 t 

OJ u 



uj 



(12) 



where h is the unit vector along the background 
magnetic field, and d e = c/oj p e is the skin depth. 

Let us now estimate the relative roles of the dif- 
ferent terms in the left-hand side of Equation [12] 
For n e = 10 11 cm -3 the electron plasma frequency 
uj pe — 10 10 s _1 , which yields the cm. The 

typical frequency of the neutral gas flows under 
consideration is uj « (1CP 2 — 10~ 3 ) s -1 and their 
length scale is L ~ 10 2 - 10 3 km. For B w 100G 
the electron gyrofrequency ujbc = 2 x 10 9 s" 1 , 
while the characteristic electron-neutral collision 
frequency v en ~ (10 s — 10 9 )sec~ 1 . Therefore the 
ratio of the second and third terms compared to 
the first can be considered small. Thus Equa- 
tion [12] can be reduced to 

t « ■ V)V£ - t(V ■ v£)}, (13) 

OJ 

which means that in the lower chromosphere the 
magnetic field is effectively frozen into the neutral- 
gas flow. 

This describes the magnetic field generation 
due to the motions of conductive fluid neglect- 
ing small scale effects scaling as the electron iner- 
tial length, and assuming the motions to be slow: 
uj << Q e . Let us now verify that the assump- 
tion made above, namely that the bulk velocity 
of the ions Vi is so close to that of the neutrals 
V n that the electric current can be written as 



~~j = n e e(vti — Ve)- It can be easily found from 
the expression [T31 that 



J 



1 

Mo 



-V x b = 



i—V x {{t -V)vt 
Mow 



(^-V)Vxt 

Mow 

To evaluate the difference of the ion and neutral 
velocities let us consider the equation of motion for 
ions: 

M ^T = e{^+[t^]}-M mn (V t -€)- (14) 
at 

As can be seen f rom Eq. [T4l the sought-after ve- 
locity deviation AVi = (Vi — V n ) is such that 
the ion-neutral drag balances the other forces ex- 
erted upon the ions. A simple estimate shows that 
the net electromagnetic force is the dominant one 
there, and Eq . [141 yields 



Ven J 



It follows then from Eq. [14] that 



MiVi, 



(15) 



and 



Al4 ~ Sv e — — < 5v e , 



where 5v e is the velocity of current carrying elec- 
trons, and ions are assumed to be demagnetized. 

On the other hand, by using Eq.[l3l the electric 
current can be estimated as 



J 



B V n 



HqL hqojHL 
The velocity of the current -carrying electrons is 
j eB V n m uj B e d\ 

Sv e —— -F7T V n < V n- 

ne fiQmujHLne^ oj HL 

(16) 

It is seen now from Eqs. [TBI and [TBI that, indeed, 
SVi «V n , and SVi < Sv e . 

A simple estimate of the currents generated at 
different altitudes and for different values of the 



S 



magnetic field can be obtained by assuming the 
neutral gas motions to be rotational. In this case 
the characteristic velocity in the vortex V n of char- 
acteristic spatial scale L rotating with the charac- 
teristic frequency lo can be evaluated as follows: 

V n ~ UJL. 

Thus, the characteristic magnetic field perturba- 
tion is 

and consequently the current density is 
SB B 



J 



fi L jioH 



One can see that the magnetic field perturba- 
tions can become quite large; the current densi- 
ties for magnetic fields of the order of 100G can 
become as great as 



Bp 



10" 



It is worth noting that this mechanism is de- 
pendent upon the velocity shear. An interesting 
feature of this process is the very high efficiency 
of the field generation around the boundaries of 
the neighbouring vortices, where the characteris- 
tic shear of the velocity of the gas motions is large. 
In this case the characteristic velocity shear is 



Vn 

L 



» cj; 



thus smaller-scale magnetic-field structures can be 
generated quite efficiently and magnetic fields gen- 
erated and current densities can become tens or 
hundreds of times higher than the estimate above. 

6. Electron and ion resistive heating. 

To evaluate the electron heating efficiency one 
should find the induced electric fields parallel and 
perpendicular to the magnetic field. This can be 
done in a straightforward way by making use of the 
component of the equation of motion for electrons 
along the magnetic field (Eq. ISJ). 



dV e \\ v en m e 

m e —r— = -eEu 

at ne 



-n 



in the lowest-order approximation one should take 
into account that the frequency u> of wave motions 



is much smaller than the electron neutral collision 
requency. This leads to the resistivity along the 
magnetic field to be classically collisional, where 
the dominant effect is due to electron-neutral col- 
lisions 



E \\ 

Bo r r 



e UJ; 



m 



.JH^Zlgt . v x {{It ■ w)t - 7f(v • %)}) 



~ v en — bB. 

To evaluate the parallel electric field one should 
take into account that the electron-neutral col- 
lision frequency varies with height from v en ~ 
5 x 10 8 to 10 10 s _1 and the plasma frequency varies 
in turn from 10 10 to 10 11 . Taking the characteris- 
tic current densities obtained above we find elec- 
tric fields in the range 

E\\ « (1CT 2 - lCT 1 ) V/m, 

with the electron heating written (Braginskii, 
1965) as 



3 , T d e T e 
2 dt 



Pe V-vt = -V-t + Qe- (17) 



Here p e is the electron pressure, and the electron 
heating Q e is determined by two "friction forces," 
one due to the relative velocity between electrons 
and ions/neutrals (R e i, R en ), and another due to 
the electron temperature gradient R T : 

with 

= -v ei m e N e {0.5{V el] - V {]l )t + (7*1 - %!)}, 

R~Z = -v en rn e N e {0^V e]r V nll )t+{v2-V n ~l)}, 
and 

Rt = -0.7N e (t,V)T e ~ |jV e ^-[7?,VT e ]. 



As a result 



Qe 



(u e i + i> en )m e N e 



{0-5^+ii+ (18) 



iV|e 2 

^(tVxit)(tv)T e }; 
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here we neglect the diference between the ion and 
neutral velocities. The last term in the equation of 
energy balance to be kept describes the heat flux 
that is written as 



<lu 



q£ = 0.7N e T e (V ell -V nll )+^N e T e -^-[t,Zl-Vil] 



and 



4.66 



We = 

N e T e {v e 



-3.16 



m e (v en + v ei ) 



-VnT fi 



2 m e il e 

We shall begin by evaluating the energy supply 
provided by resistive dissipation that is deter- 
mined by first two terms in the right hand side 
of the equation of energy balance. The volumetric 
heating power of the electrons can be estimated as 



3 t dT - ■ w 
~nk B —~ nEll 



Mo 



L 2 



^(io- 4 ^io- 3 )-^. 

In a thermal equilibrium when the macroscopic 
flow evacuates the power supplied by collisional 
resistive dissipation, one can estimate the temper- 
ature variation with altitude as 



dt 



V, 



die die 

dz ' dz 



\nk B V e: 



where V ez is the characteristic macroscopic verti- 
cal velocity of electron (plasma) motion that can 
be evaluated to be of the order of the sound veloc- 
ity. Then 



die 

dz 



j\\ E \\ 
lnk B V s 



(0.01 - 0.1) V/kr 



It is easy to see that the collisional heating of 
the electrons becomes rather efficient. Another 
possible evaluation of the temperature variation 
with altitude can be obtained by assuming a sta- 
tionary static thermal equilibrium. Under such 
conditions the power supply will be balanced by 
the divergence of the heat flux in Equation IT7l 
which can be estimated as 



77> 



dz Kz dz 



where k z is the electron thermal conductivity in 
the parallel direction. For our conditions when 
electron-neutral conditions are dominant in the 
parallel direction, 



n P T R 



*z ~ 7 , 

where 7 is a coefficient of order unity (according to 
Braginskii, 1965) for electron-ion collisions; when 
ions are singly ionized it is equal to ^3.8. This 
leads to the following estimate of the characteristic 
heat loss rate 



_v d k B n e T e dk B T e 

ql « ^-7 5 — 

oz m e v en oz 



k%n e T 2 

T 2 1 



where (T e /Lx) is the characteristic variation of 
temperature with altitude. Under the physical 
conditions of the chromosphere this will lead to 



, T e ^ j\\E\\m e v en 
L T ki,n e 



1/2 



0.1 - 0.3 eV/km. 



Taking into account an exponential decrease 
of the plasma density one can find that it can 
heat plasma to the hydrogen ionization energy 
(13.6 eV) on a characteristic distance of the order 
of a few hundred km. It should be noted however 
that the thermal conductivity increases with the 
temperature and this can result in some flattening 
of the temperature profile. 

Another important physical effect is related to 
the ion heating due to the perpendicular electric 
fields. It is known that the major component of 
the current along the electric field perpendicular 
to the magnetic (the Pedersen current) can be car- 
ried by ions. This can lead to the more efficient 
heating of ions than of electrons. To evaluate the 
Pedersen current we must estimate the ion veloc- 
ity deviation from the velocity of neutrals in the 
direction of the electric field, as was already done 
in the previous paragraph 



A 



3]> 



Mm, 



1 



MiV m n 
It follows then that 



[3 x 



2*1 



3 Bp 

nMiVi. 
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nMiV m m e v in L 2 

Ion heating due to the perpendicular electric held 
component can be evaluated as 

— l - ~ (JZ ■ {£ + [V n x t ]}) ~ -i — s--^ 

This estimate shows that the ion collisional 
heating can be several times smaller than, or com- 
parable to, the electron heating. 

7. Discussion 

The current generation mechanism that we pro- 
pose here is strongly dependent upon four basic 
parameters: the velocity of the turbulent gas mo- 
tions, the background magnetic field, the char- 
acteristic frequency of the spectrum of turbulent 
motions of the gas, and the characteristic spa- 
tial scale of velocity shear. The physical process 
of the generation of electric currents and mag- 
netic fields is actually nothing other than the well- 
known turbulent-dynamo mechanism operating as 
a result of the the freezing of the magnetic field 
into the flow of the neutral gas. 

The altitudes where this mechanism will oper- 
ate are strongly dependent upon the local mag- 
netic field strength: the stronger the field, the 
lower the altitude of electron magnetization. The 
efficiency of the electric current generation is pro- 
portional to the magnitude of the velocity shear. 
The characteristic spatial scales are proportional 
to the characteristic scales of this velocity shear, 
i.e. the flow vorticity. It is interesting to note that 
in this formalism, magnetic fields and currents 
can be generated by both compressible and incom- 
pressible motions (see the first and second terms 
in equation [T"3]) . Compressive motion tends to 
amplify or weaken the background magnetic field, 
while the rotational component effectively gener- 
ates the perpendicular or helical components of 
the field. This means that compressional motions 
tend to generate transverse currents, while non- 
compressional motions can generate field aligned 
currents as well. 

Since density tends to decrease exponentially 
with height in the chromosphere, the character- 
istic velocity of chromospheric flows can increase 
by a factor of ten or more over their photospheric 



values. This gives rise to an estimate of the charac- 
teristic magnetic field generated by such motions 
readily becoming comparable to or even larger 
than the background field. Their relative mag- 
nitude is determined by the ratio of the charac- 
teristic velocity of neutral motions at the altitude 
of electron magnetization to the characteristic pa- 
rameter luH. One can see that for characteris- 
tic frequencies of the order of 10~ 2 — 10~ 3 s _1 , 
and characteristic heights of the order of 140 km, 
the field generation becomes quite efficient already 
for the velocities of the order of ten-several tens 
of km/s. The helical magnetic field component can 
become larger than the background field. A sim- 
ilar estimate shows that amplification/ weakening 
of the background magnetic field component can 
be evaluated to be of the order of 

b Sn 
B n ' 

""Under such conditions, we show that the elec- 
tron heating can become rather efficient due to 
the collisional resistive dissipation of the electric 
current. These processes take place in the regions 
where the magnetic field vorticity increases. 

It is worth noting that the efficiency of the 
electron and ion heating can become substantially 
more important than that due to the collisional 
resistive dissipation described above. The electric 
current and magnetic field generation we have dis- 
cussed takes place at the chromospheric footpoints 
of flux tubes that may extend into the corona. The 
direction of the generated field is determined by 
the characteristics of the gas velocity shear and 
can lead to the formation of the magnetic field con- 
figurations that can become unstable, or can form 
new structures by interaction with neighbouring 
fields with different orientations. The increase of 
the perpendicular component of the magnetic field 
can result in kink instability if this component sat- 
isfies the Kruskal-Shafranov condtion 

SB I 
~B ' > R' 

Here I is the characteristic scale of the cross sec- 
tion of the tube and R is its characteristic length. 
Under such conditions the tube cannot keep its 
configuration intact; it will become unstable and 
then different kinds of perturbations will result in 
its reconfiguration. Other types of emerging con- 
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figurations with magnetic field inversions will re- 
sult in local reconnection. The detailed study of 
these effects lies beyond the scope of our paper 
and will be carried out elsewhere. 

The result of the heating would be a rapid 
increase of the degree of ionization of the gas. 
It is worth remembering that our whole descrip- 
tion is valid only when the degree of ionization is 
low. This rapid ionization process can lead, ac- 
cording to most radiativ e-convective models (e.g., 
Stein fc Nordlundlll998h to an efficient decrease of 
the radiative cooling rate of the gas. Thus the gen- 
eration of electric currents results not only in the 
plasma heating itself, but also causes a decrease 
of the cooling process that will then strengthen 
the effect of radiative heating. We cannot make 
any quantitative estimate of these effects, but they 
may be adressed in future simulations. 

The coupling of acoustic-type waves and MHD 
waves can also be treated in th e MHD approxima- 
tion fe.g.. iHollweg et al.lll982i ). The major differ- 



ence in our approach here is that the neutral colli- 
sions reduces the feedback effect on the motion of 
the gas resulting from the newly generated mag- 
netic field. This produces a more efficient transfer 
of energy from gas motions to current generation 
and magnetic field amplification. 

It seems possible to incorporate certain aspects 
of this formalism into a system of conservatio n 
equations similar to those used in lAbbettl (|2007| ). 
In that work, the resistive MHD equations were 
solved numerically within a computational do- 
main that includes both a turbulent model convec- 
tion zone and corona. We are currently updating 
this system of equations to incorporate additional 
physics, and are planning to perform a compar- 
ison between a standard resistive radiative-MHD 
model and our new results. We hope to report on 
the results of this study in the near future. 

8. Conclusion 

We have presented an analysis of the effect of 
electric current and associated magnetic fields gen- 
eration in a "chromospheric dynamo layer" where 
electrons become magnetized while ions remain 
collisionally coupled to the neutrals. We have 
shown that electric currents and magnetic fields 
can be generated very efficiently due to turbulent 
motions of the neutral gas. The efficiency of this 



physical process is proportional to the character- 
istic velocity shear of the gas turbulent motions 
and inversely proportional to their characteristic 
frequency. We also found that the magnetic fields 
thus generated can be comparable to or even larger 
than the background magnetic field for motions 
having characteristic scales of the order of several 
hundred km and characteristic time scales of the 
order of several minutes. This can produce a sub- 
stantial restructuring of magnetic field configura- 
tions and an opportunity to create multiple sites 
of reconnection. It also results in an efficient in- 
crease of collisional resistive heating of electrons 
and ions, and hence the rapid ionization of the gas, 
thus unstably altering its thermal equilibrium. 
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